Ashkin- Teller model on the iso-radial graphs 

S. Lotting M. A. Rajabpour 2t 

1 Westfalische Wilhelms-Universitat Miinster, 
Wilhelm-Klemm-Strafie 9, 48149 Miinster, Germany 
2 Dip. di Fisica Teorica & INFN, Universita di Torino, 
Via P. Giuria 1, 10125 Torino, Italy 

April 1, 2010 
Abstract 

Wc find the critical surface of the Ashkin- Teller model on the generic iso-radial graphs 
by using the results for the anisotropic Ashkin-teller model on the square lattice. Different 
geometrical aspects of this critical surface are discussed, especially their connection to the 
anisotropy angle. The free energy of the model on the generic iso-radial graph is extracted 
using the inversion identities. In addition, lattice holomorphic variables arc discussed at 
some particular points of the critical line. We check our conjectures numerically for the 
anisotropic triangular-lattice Ashkin- Teller model. 
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1 Introduction 

The Ashkin- Teller model, as a simple spin model exhibiting a critical line in its phase diagram, 
has been studied for many years. This model was introduced by Ashkin and Teller [I], and 
can be seen as two coupled Ising models. Some physical realisations of the Ashkin- Teller 
model are layers of molecules adsorbed on surfaces, such as selenium on nichel, molecular 
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oxygen on graphite, and atomic oxygen on tungsten [21 [31 0]. The model can be mapped, 
on the square lattice, to the staggered 6-vertex model, and the phase diagram and some of 
the critical exponents were found by mapping the model to the Solid-on-Solid (SOS) model 
[3 El [3 E] • Although the phase diagram is qualitatively quite well known, the model has not 
been solved exactly for the generic coupling constants as the staggered 6-vertex model has 
not been solved so far. The staggering of the model disappears on the critical line, where 
one can completely solve the model [7]. The quantum version of the Ashkin- Teller model 
was also studied by considering the highly anisotropic model on the square lattice and has 
a different phase diagram [9]. The continuum version of the model is as much complicated 
as the discrete version. It is related to the sine-Gordon model [9] on the critical line, see 
[ID] for more recent developments. On the critical line one can also write the action with 
respect to fermionic variables (it is the massless Thirring model [H]). The corresponding 
conformal field theory of the model on the critical line, called orbifold conformal field theory, 
was discussed in [12} [T3l W$\ [T5l [T6] , see also [T7] . The four-point correlation functions of the 
spin operators were also calculated exactly by conformal field theory methods [H] . Although 
in the physics community it is widely accepted that the above continuum models describe 
the critical line of the Ashkin- Teller model, there is no rigorous proof for the validity of the 
results. 

There was no systematic way to study the continuum limit, even for simpler systems 
such as percolation or Ising models, before the invention of Schramm-Loewner evolution by 
Schramm [19], for review see [20]. This opened the way to investigate the continuum versions 
of the critical lines in statistical models. Recently Chelkak and Smirnov showed that this 
method is also useful to investigate the universality of statistical models. They proved that 
the Ising model on different iso-radial graphs, which will be discussed in Section [3l have 
the same continuum limit at criticality [21]. The iso-radial graphs are important because 
they are the most generic graphs preserving some properties of complex analysis and have a 
well-defined continuum limit. They are also the most general graphs where one can expect 
conformal invariance, in the thermodynamic limit, for critical statistical models. 

These reasons are enough to stimulate theoretical physicists to investigate statistical mod- 
els on iso-radial graphs. The study of statistical models on planar graphs was initiated by 
Baxter [22\ [23] . He introduced the most generic planar graphs such that the star-triangle re- 
lations (or Yang-Baxter equations) are well-defined for the eight vertex model. These graphs, 
now called Baxter lattices, are those made of the union of simple non self- intersecting curves 
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crossing the plane, with the property that a given curve can intersect any of the others at 
most once. Iso-radial graphs are a special embedding of the Baxter lattices. Baxter solved the 
eight- vertex model [22] and the Potts model [23] on these graphs. With the above motivations 
about the importance of iso-radial graphs, we will study the critical line of the Ashkin- Teller 
model on them. Most of the studies on the Ashkin- Teller model are performed on the square 
lattice. To the best of our knowledge the formula for the critical line of the Ashkin- Teller 
model on the generic iso-radial graph is not written anywhere. 

In this paper, by using the formula for the critical line of the anisotropic Ashkin- Teller 
model on the square lattice, we will find the critical surface on iso-radial graphs. The central 
claim is, the formula for the critical point of the model on the anisotropic square lattice is 
sufficient to find the critical point on the most generic iso-radial graph. This was already 
checked for the Ising model in [24] and for the Fateev-Zamolodchikov point of Zjv models [25] 
in [26]. Here we will investigate its validity for the critical line of the Ashkin- Teller model. 
The paper is organised as follows: in Section [2] we define the Ashkin- Teller model and recall 
its different phases on the square lattice; in Section [3] we first review the general method for 
the Ising model with arbitrary couplings on iso-radial graphs; in Section [H by employing the 
method of Section [31 we introduce our conjecture for the critical line on iso-radial graphs, 
then, in Section [5l we provide numerical calculations for the anisotropic triangular lattice 
in support of it. Finally, in the last Section we will conclude our findings and also discuss 
possible future investigations. 

2 Ashkin- Teller model: definition 

The partition function of the Ashkin- Teller model is given by: 

* = En^)=En^> 

C i,j C i,j 

Sij = fiipiOj +TiTj) + a{pi(JjTiTj) , (2.1) 

where <7j and T{ are Ising variables, i and j run over all pairs of neighbouring sites, and the 
summation is over all possible configurations. One can write the above partition function 
with respect to complex variables as follows: 

z = e n w &i) = e n coo- + xis ^- + x ^ s2 j + xi ^ 3 4) . ( 2 - 2 ) 

C i,j C i,j 
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where Sj = ai ^' and 



Co = 1 + tanh a tanh 2 (3 , (2-3) 

e 2/3 _ e -2/3 

Xl = e 2/3 + e -2/3 + 2e -2a ' ( 2 - 4 ) 

e 2/3 + e -2 /3 _ 2e -2a 

X2 " e 2/3 + e -2/3 + 2e -2a • W 

In this notation the Z4 symmetry of the Ashkin- Teller model is evident. Since some of the 
forthcoming formulas have a simpler form in one or the other of the above notations, we will 
mostly provide them in both notations. The duality transformation in the /3, a notation is 
implemented by the relations: 

1 / 1 + tanh 2 /3 tanh a - 



~ . i / 1 + tann a tann a \ 

pip, a) = -log g 

y 1 4 6 V tanh 2 /? + tanha / 



tanh j3 + tanh a 
1 / (coth /3 + tanh /3 tanh a) (coth j5 + tanh ft coth a) . 
° (/5 ' a) = 4 l0g ( 2 + tanh a + coth a ) 

In the x\, X2 notation they become simply: 

1 — x 2 



Xl 



I + X2 + 2xi 



1 + x 2 - 2xi /„ 7A 

x 2 = — — — • (2.7) 

1 + x 2 + 2xi 

On the arbitrary lattice, the physical region of the model (i. e. that with positive Boltzmann 
weights) is the triangle enclosed by the lines 

1 + 2xi + x 2 > , 1 - 2xi + x 2 > , x 2 < 1 . (2.8) 



To stick to the ferromagnetic region we also need to consider x\ > 0. From Eq. (I2.7P one can 

easily find the self-dual line xf + 2xf = 1; on the square lattice, this line is critical for x\ > 3. 

On the square lattice one can write the following equations for the correlation exponent [7]: 

12 , . 

- = 2--, (2.9) 

^(x) = (210) 

It is useful to identify some well-known points on the critical line. The first interesting point 
is the 4-state Potts model. On the square lattice it is given by the following values in the 
(3-a and x\-x 2 planes: 

^ log 3 ~ 0.274653 , (2.11) 

(2.13) 
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The next interesting point is the Fateev-Zamolodchikov (F-Z) point [25] . with the following 
values: 



f FZ = -I log f 17 - 8\/2 + (8 - 2V2){J 2-V2-J2 + V2)) ~ 0.302923 



a£ z = - log(\/2 + 1) ~ 0.220343 , 



*ifz = — t^T - 0.351153 , (2.14) 



sin (lg) 
sin(ff) 

« sin(^) sin(ff) 

4fz = 7^7 — t^- 0.297693, 

sin(ff)sin(If) 

3 

VFZ = V 

The third interesting point is the Ising point; it is the critical point along the Ising line a = 
(or X2 = xf), and has the following values: 

Pf = \ log(l + y/2) ~ 0.440687 , of = , 

xf T = V2-1, x$j = 2-V2, ( 2 - 15 ) 
i// = 1. 

The next interesting point is the one lying on the boundary of the physical region; we call it 
the "terminator" point: 

$ = + log V2 -> 00 , (2.16) 

xf T = i, xf T = 0, (2.17) 

z/ T = 2 . (2.18) 

On the square lattice, the X-Y point (where v — > 00) is not in the physical region: if we 
naively demand this value of the correlation exponent from the equation (I2.9D then we will 
get: 

* 1 c 2sin(f)-l 

x° XY = — ~ 0.566454 , xf xy = ^ ~ -0.132909 , (2.19 

1XY l + 2sin(f) ' 2XY 2sin(f) + l ' v ; 

but the translation to the /3-a plane poses some problems, as one gets complex-valued cou- 
plings: 
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P S XY = \ log [- (1 + 2^2 + 2(1 + ^)^2- V2)} , a s XY = ilog(-(V2-l)) . (2.20) 
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We will conclude here the introduction of particular points on the square-lattice critical line; 
in Appendix [A] we will locate those points on the triangular and honeycomb lattices and 
provide more details. 
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Figure 1: The phase diagram of the Ashkin- Teller model. In region I we have (s) 7^ and 
(s 2 ) 7^ 0, in region II we have (s) = (s 2 ) = 0, and finally in region III (s) = and (s 2 ) 7^ 0. 
The blue lines (mapped onto each other by duality) are Ising lines separating region II from 
I and III. Thick black lines are the boundary of the physical region and the red line is the 
self-dual line: it is critical between points P and T. I Q , I\ and I2 are Ising points. 
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3 Ising model on iso-radial graphs 



In this Section we review the method employed in \24\ I26j to extract the critical point of the 
Ising model on iso-radial graphs starting from the critical point of the anisotropic square- 
lattice Ising model. Here we generalize the situation to the off-critical region as well. Then, 
in the next Section, by using the solution for the anisotropic square-lattice Ashkin- Teller 
model, we will derive the critical surface of this model on the generic iso-radial graph. 

To see the importance of the anisotropy angle in the Ising model, consider a square lattice 
with different couplings in the x and y directions. To fix the notation write the partition 
function as: 

Z = ^2 Wi(x,y,x',y') = i 1 + w xS x , y s x ', y + w y s X) yS x>yl ) , (3.1) 

C x,y,x',y' C x,y,x',y' 

where the primed coordinates denote nearest-neighbours of (x,y). In this case, in order to 
maintain invariance under lattice translations and reflection symmetry about the x and y 
axes, the only allowed transformations are relative scalings of the x and y [27]. In our case 
each elementary square of the covering lattice is deformed into a rectangle. To proceed, we 
define disorder variables fj,? on the sites of the dual lattice, which is another rectangular lattice 
with vertices at the centers r of the faces of the original lattice. The insertion of a disorder 
variable corresponds to modifying the weights so that the order variable s r has monodromy 
s r — s r on taking the point r in a closed circuit around f. This is equivalent to the 
introduction of a path on the sites of the dual lattice from f to infinity, such that the weights 
on edges (rr') intersected by the path are modified by the substitution s r s r > — > —s r s r >. To fix 
the gauge we will consider a path parallel to the x axis and going to +oo. Thus the disorder 
operator has the following form for the Ising model: 

, f = n [?; Wr i 0.2) 

x x (1 + W rr 'S r S r i) 

(rr'Jintcrscctcd by path 

It is easy to see that the following equation is valid for two adjacent disorder variables: 

1 - w y sis 2 

Mi = t- H2 • 3.3) 

1 + W y SlS2 

Now consider the quadrilateral whose vertices are the two neighbouring spin variables s± and 
S2 and the two neighbouring disorder variables /Uj and /xs: this quadrilateral is living on the 
covering lattice. Using Eq. (|3.3p one can show that: 

- + e-^- 6 H 2 i + V 2 2 + e% S = + ^T^ 1 ) ^ ~ ifci) , (3-4) 



Figure 2: The faces of the covering lattice form rombi with angle 6; the thick dashed line is 
the frustration line. 

with spinor-like variables tp defined as: 

■011 = e l §si^j , V 2 i = e-**82Hi > ,„ rN 

■ (T -fl) ( 3 - 5 ) 

^12 = e * 2 *lA*2 > ^22 = S 2/i2 • 

To give a geometrical meaning to 9 consider the quadrilateral (si, //j, S2, to be a 
rhombus with angle 6 on the spin-variables vertices. Then one can write Eq. (13. 4p as follows: 

^2^ e 5z e = [e %e H ^'^l — ijj 2 i) , with ij> r ? = e ^ s r ^,f , (3.6) 

where is the angle that the directed segment rf makes with the axis parallel to Si-//j. The 
left-hand side of the above equation is akin to a discrete contour integral around elementary 
rhombi of the covering lattice. By replacing % — >■ —i one can find a similar equation for the 
discrete antiholomorphic spinor-like variables, with the following property: 

-f + ^e ^ 2 y - 

J y 1 + <V 

Using the Eqs. (|3.6p and (|3.7p one obtains the following equation: 



#11 + + *2i + <fe = (> + + ^Tfy <*21 - «*12> • (3-7) 



S> fe « = < 2sta ©rr5 - 2cos (I)if^j) E* • < 3 - 8 > 

e » we 

The above relation is clearly similar to the equations of motion for a field theory of a free 
Majorana fermion, i. e. : 

r t / 

dzX = , (3.9) 



S 



Figure 3: 2-colorable Baxter lattice and its Kenyon and Schlenker rombus embedding. We 
define the Ising variables in the shaded region and the disorder variables in the white region. 



where x an d X are the chiral components of the free Majorana fermion. The comparison then 
yields: 

X = aip , x = H , a = ib , (3.10) 

/0\l-wl /0\ 2w v 
m = 8sin(- |-8cos(- ^. 3.11 

It is easy to see that m = for w y = tan(|), and that tp is the discrete holomorphic operator. 
This means that this particular embedding of the covering lattice leads to a critical system. 
This gives a complete geometrical meaning to Eq. (|3.9p . The discrete Dirac equation for the 
Ising model was derived long ago by many authors, see for instance [28J . 

It is possible to extend the above results to the more general graphs called iso-radial 
graphs. To introduce iso-radial graphs we need first to define the Z-invariant or Baxter 
lattice. This is a planar graph C which is a union of M simple (non-self-intersecting) curves 
crossing the complex plane from Xj — too to x'j + ioo, where the values Xj and x'j (1 < j < M) 
are distinct, and with the further property that a given curve can intersect any of the others 
at most once (Fig. [3]). The faces of C are 2-colorable, and in general we require that its 
vertices are all of degree four. Consider now the planar graph Q whose vertices are associated 
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to shaded faces of C. We can define an Ising model on Q with general weig hts xf. The 



vertices of the dual lattice Q* then correspond to the white faces of C The vertices of the 
covering graph C (the union of the vertices of Q and Q*) correspond each to a face of C, 
irrespective of color. There is a theorem, due to Kenyon and Schlenkcr [29J, stating that 
there is a rhombic embedding of C into the plane, that is, one in which all the edges have 
equal length, see Fig. [3j Each rhombus corresponds to an edge E of G, and defines an opening 
angle Be- This graph is called iso-radial graph because one can draw circles with equal radius 
centered on vertices of Q; the intersection between the circles will then define Q* , see Fig. [5J 
These kinds of graphs are very important because they are the most general graphs that can 
offer a discrete complex analysis with a well-defined continuum limit. Most of the well-known 
planar graphs are just instances of iso-radial graphs. For example, the triangular lattice is 
an iso-radial graph with 9 = ^ and the honeycomb lattice has 6 = As before, we can 
define the Ising model and its dual on this graph and extract Eq. (|3.8p with the substitution 
6 — > 6e- It is easy to see that one can also extract the critical point of the Ising model on 
iso-radial graphs by requiring: 



The above simple example shows that one can extract the critical line of some models 
starting from the sole knowledge of the critical point for the anisotropic square lattice. In the 
next Section we will use the same technique to extract the critical line of the Ashkin- Teller 
model on the general iso-radial graph. The above argument shows that the the critical and 
off-critical Ising model may have a well-defined continuum limit on iso-radial graphs. The 
continuum limit of the critical Ising model was already rigorously investigated in |21j . and 
the above argument indicates that possibly the same graphs give the right answer in the 
continuum limit for the off-critical case. 

4 Ashkin- Teller model on the iso-radial graphs 

In this Section we will first give the critical surface for the Ashkin- Teller model on the generic 
iso-radial graph and then we will study holomorphic variables of the model at some particular 
points. 

4.1 Critical line of the Ashkin- Teller model 

In this Section, using the same argument as in the previous one, we give the exact form of the 
(integrable part of the) critical surface of the Ashkin- Teller model on the iso-radial graphs. 




(3.12) 
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Figure 4: The corresponding iso-radial packing of the Kenyon and Schlenker rombus embed- 
ding. At the critical point Be is the only parameter needed to describe the critical temperature 
of the system. 



The only formula we need for this purpose is the critical line of the anisotropic Ashkin- Teller 
model on the square lattice. It is given in |30j and has the following expression: 

= (4-1) 
tanh^AJ , (4.2) 

tanh(2/3j / ) , (4.3) 



sin(u) 
sin(A) 
sin(A — u) 
sin(A) 



A sinh(2a 2 ,.) sinh(2aj,) 

- = cos(A) = sinh(2Ac) = sinh ( 2j9y ) ' ^ A > 

where B is the anisotropy angle, u is the spectral parameter and A is a parameter labelling 
the universality class^. The above formulas give the two-dimensional exact solution manifold 
of the anisotropic Ashkin- Teller model among the four couplings. By "exact solution" we 
mean that they satisfy the star-triangle relations on this manifold [54] . This solution can 
not be the full critical surface of the Ashkin- Teller model, however, it is the largest manifold 
satisfying the star-triangle relation. The parameter A has the following lattice-independent 



1 Specifically, we have Xp = 0, Xfz — f , \i — f , At = §7T and \xy = §7r. 
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relation with the correlation exponent (as can be obtained from the square-lattice case): 

2vr-2A 

Now the only thing we need to do is considering 9 as the angle of the rhombus in the 
iso-radial graph. Let's consider some examples: for the square lattice we have = 5, which 
leads to the following critical line: 

sinh(2/3 s )e 2QS = 1, (4.6) 



x 



1 - 2xf , (4.7) 



and, for the parameter A: 



A 9 = 2 arccos — - . (4.8) 

ft 2tanh2/3 v ; 

To get the critical line of the triangular lattice we need to consider 9 = ^; then we have 
the following lines in the f3-a and x\-X2 planes: 



= 4a T (p ip T 



1) = 2, (4.9) 



X T = 1 + X T _ ^ 5 ( X T )2 + 4x T ( 41Q) 

Here A has the following relation with the coupling constants: 



/ /3sinh2/3-cosh2/3\ 

At = 3 arcsin * / — . (4.11) 

I y 4sinh2/3 J v ; 

It is easy to see that in this case one can reach the X-Y point v — > oo by staying in the 

physical region, i. e. with real f3 and a. A similar calculation can be done for honeycomb 
lattice, with 6 = t?, finding: 

e 4Q %inh 2 (2/? H ) - e 2 "" cosh(2/3 H ) = 1, (4.12) 

xf = 1 - 2(acf f . (4.13) 

In this case A can be written as: 



/ cosh(2/3) + VScosh 2 2/3 -4\ 

A/f = 3 arccos , , — ^ • (4.14) 

I 4sinh2/3 i v ; 

On the honeycomb lattice, the universality classes falling in the physical region are a subset 
of those accessible in the triangular and square lattices (for a comparision of the three cases, 
see Fig. [5]). The same conclusions can be extracted by using the duality relations, Eqs. (|2.6p 
and (12. 7J) . The above results are consistent with those in 1311 1 2 ! 



2 Note that some parts of that paper, identifying the critical line of the Ashkin- Teller model with the 
universality class of the Potts model, are wrong. 
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A-T phase space, fS— a 




-0.8 I 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

A-T phase space, Xj— x> 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



X l 

Figure 5: The critical lines for the Ashkin- Teller model on the triangular, square and honey- 
comb lattices (denoted by the symbols A, □ and O respectively), depicted in the /3-a (top) 
and the xi-x% (bottom) planes. Dashed lines are the Potts and Ising lines; the critical lines 
span the interval between the Potts (P) and the X-Y (XY) points; the self-dual line is also 
shown outside the (square-lattice) critical interval. In the %\-xi plot, the shaded region is 
the non-physical one, i. e. the one having negative Boltzmann weights: as such, that region 
cannot be mapped to the /3-a plane, and there the points on the threshold are mapped to 
infinity (they are the X-Y point for triangular, the "terminator" T point for square and a 
point labelled U for the honeycomb lattice). For more information on the triangular- and 
honeycomb-lattice critical lines, see Appendix [XI 
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The same procedure is also effective for the anisotropic triangular lattice with different 
weights f3i and ai on the triangle sides i = 1,2,3. Consider a (proper) embedding of the 
regular triangular lattice into the plane by some general linear transformation of the coor- 
dinates. Then locate the dual vertex at the circumcenter, that is, the point at which the 
three perpendicular bisectors of the edges meet, equidistant from the three vertices. This 
construction guarantees that adjacent pairs of vertices always form a rhombus. Each triangle 
is associated to three different rhombi of angles 0\, 62 and #3 = tt — 6\ — 62- Finally, one can 
write the critical line as follows: 



sin (A — Ui 
sin (A) 



" (A = ^ , (4.15) 
A 

tanh(2/3i) , (4.16) 



A . sinh(2aj) , A „„. 

T = cosA = g - ,L > 4 - 17 

2 smh(2pij 
with i = 1,2,3. It is not difficult to see that the above formulas give a three-dimensional 
manifold among the six thermodynamic parameters. This can not be the full critical manifold 
of the system. The point is, our method of extracting the critical manifold is in close relation 
with the integrability of the system. On this manifold the star-triangle relations among the 
couplings are satisfied, meaning that our method works only as long as the star-triangle 
relations are satisfied; the same observation was also made for the Fateev-Zamolodchikov 
point in the lattice Z^v models [26]. The close connection of our method to the star-triangle 
relations can be related to the integrability properties of the iso-radial graphs and can come 
from deep geometrical reasons, see |32] and references therein. The above formulas are also 
true for the anisotropic honeycomb lattice. The only subtlety is, we need to consider three 
angles 9\, 62 and #3 constrained by 6\ + 62 + #3 = 2-tt. 

4.2 Free energy 

In this Subsection we write the exact free energy of the Ashkin- Teller model on the critical 
surface of iso-radial graphs. To avoid any non-positive Boltzman weight we will concentrate 
on the region between the Ising and the 4-state Potts models. To find the free energy on 
iso-radial graphs we first need to find the free energy on the anisotropic square lattice; in 
order to do so, we need the unitarity conditions [33]. We follow the notation of [30] to keep 
the equations simple. Consider the weights 

W x (i,j) = p x exp[f3 x (aiaj + TiTj) + axiaiVjTiTj)} ; (4.18) 
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W y (i,j) = p y ex.p[f3y(aicrj + TiTj) + ayiaicrjTiTj)} , (4.19) 

where p x and p y are some normalisation constants, that can take any finite value and only 
shift the free energy by a constant. It is convenient to use the following parametrisation: 

_ sin(ix) _ sin(A - u) _ sin(A + u) 

sm(A) sin (A J sin (A J 

If we choose p X:V = -^ eax ' v tanh 2/3 XjV , the unitarity conditions, or inversion identities, will 
be: 

Y,W x (i,k\u)W x (k,j\-u) = 2s+s-6(i,j) , (4.21) 
k 

W y (i,j\u)W y (i,j\-u) = 2 S + S -- ( 422 ) 
Following \22\ [33] , one can argue that the free energy should satisfy: 

f{9) + f(n + 9) = log(2 S+S _) , f(9) = /(vr - 9) . (4.23) 



The above equations were solved in [34]; one can write the free energy of the Ashkin- Teller 
model on the anisotropic square lattice as follows: 

/(M) = log( SS _) + 2 r Slnh[( r 7 tanh(Ax) cosh \( 1 - 2-)\x] — . (4.24) 
Jq sinh(7rx) L\ it J J x 

For the free energy of the system on the iso-radial graphs consider a graph with a very 
large number of bulk edges Ne, so that the number of exterior sites scales as \/Ne~: then, 
by following Section 5 of [22], one can show that the free energy for the iso-radial graphs is 
simply the sum of the per-site free energies: in other words, 

/= E f(0ij) + O(y/N^), (4.25) 

fo'GS(0)) 

where the function f(0ij) is the same as in Eq. (|4.24p . and the summation runs over all edges 
of the graph Q. It is worth mentioning here that to get the free energy of the system on 
iso-radial graphs it is sufficient to know the free energy of the system on a particular lattice, 
not necessarily square; in other words, /(%) does not depend on the specific details of the 
lattice. 



3 The free energy of iso-radial graphs is the same as in Baxter lattices. 
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4.3 Holomorphic variables at some particular points 

In this Section we want to show that it is possible to extract some holomorphic variables, 
written with respect to spin and disorder variables, for the Ashkin- Teller model at some 
particular points on the critical line. The holomorphic variables can be written either in 
terms of a and r or s variables. The definitions of disorder variables with respect to the a 
and r spin variables are: 

l + x 1 {-a 1 a 2 + TxT 2 ) . . 

^1 = 1 j- Z („ ~ X ~ ~ r Ho5 > ( 4 - 26 ) 

l + ZiOlCr 2 -TiT 2 ) , , 

Mrl = T~ 7 1 ■ ( 4 - 27 

1 + xi(aia 2 + t\Ti) 



Vsl = -i , Z t„ T* i JTi , I /„*2^ ^2 , (4.28) 



The disorder variable with respect to the complex spin variable s is given by: 

1 + Xi(uj*SiS* 2 + WSJS 2 ) - CJ 2 X 2 (s* 2 S 2 ) 

1 + xi(sis* 2 + s*s 2 ) + x 2 (sfsl) 
where w 4 = 1. There is no simple formula directly relating [l s and (fi a ,fi T ). From the first 
two equations one can get, for the Ising point, the following two parafermionic operators on 
the generic iso-radial graphs: 

VVrf = e 2~a r fj> a f , tprrf = e ^T r p, T f , (4.29) 

where, again, 9 r ? is the angle the directed segment rr makes with the axes parallel to ai-fj, a f 
or Ti-/i|. It is easy also to show that at the 4-state Potts point the following variables satisfy 
the discrete Cauchy-Riemann equations on the iso-radial graphs: 

The above discrete variables could be candidates for the continuum counterparts with spin 
p = 1. Using the variables s and fi s2 one can show that the following quantity satisfies, at 
the F-Z point, the discrete Cauchy-Riemann equations on the iso-radial graph: 



4>srr = e ^ sffif . (4.31) 

The above holomorphic variable is the discrete version of the parafermionic operator of 
Ref. [25]. The next interesting holomorphic variable is related to the X-Y point and has 
the following form: 

ip srf = e~ lJ Tr sin? ■ (4.32) 

The above variable can be the discrete version of the parafermionic operator of spin p = \ 
in the continuum limit. From the continuum limit of the critical line of the Ashkin- Teller 
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model, we know there should be some holomorphic operators all the way along the critical 
line [6]. The spin of these holomorphic operators is given by p = |, with r given by Eq. (|2.9p . 
The above four operators are consistent with this continuum prediction; however, it seems 
that at points other than those four it is not easy to write the holomorphic variables as a 
function of the discrete spin variables. 



5 Numerical results for the anisotropic triangular lattice 

This Section presents Monte Carlo results in support of the formula for the critical manifold, 
Eqs. I4.15ll4.171 in the particular case of a triangular lattice with couplings f3±, a± on two sides 
of the triangle and /?2,«2 on the other. Note that, in the final formulation of the model (as 
well as in the simulation algorithm used), only the values of the couplings on the different 
links encode the geometrical anisotropy and carry the information about the opening angle 
of the rhombi. 

On such a system the angles are constrained by 28 1 + 62 = vr; in the four-dimensional 
coupling space, the full critical line is expected to be three-dimensional, while its integrable 
part (in the sense specified earlier) is a two-dimensional surface. After some manipulation 
of the seven starting relations (the constraint on the angles, plus Eqs. I4.15H4.171 taken for 
i = 1 and i = 2), one can express the critical surface £ as explicitly parametrised by the two 
couplings /?i , a\ : 



1 /l - sinh z 2/3i + 2sinh z 2ai - 2sinh2ai cosh2ai 

Bo = -arctanh ^ 

H 2 V cosh 2 20i 

= G03i,ai); (5.33) 

1 /sinh 2G({3i, a\) sinh 2a\ \ 
-arcsmh — 

2 V sinh2/3i / 

= H(p liai ). (5.34) 



This surface is the union of one-dimensional "iso-A" lines, where the universality class stays 
constant; the universality class label is still given by Eq. 14.171 

In order to confirm the validity of the above formulas, we chose four highly anisotropic 
points in the critical region between the four-state Potts and the Ising "iso-A" lines and 
sampled the (density of) energy susceptibility along segments crossing the critical manifold. 
Then, the scaling behaviour of this quantity can be checked against the predicted from 
Eqs. 14.171 and 14.51 Such a signal is maximised if the segments are chosen to lie orthogonal 
to the critical surface, whose orientation can be found analytically from the parametrisation 
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Point name 


Pi 




ft 


a 2 


expected v 


Al 


0.205000 


0.202493 


0.115419 


0.113956 


0.678871 


AFZ 


0.250000 


0.180302 


0.085211 


0.060398 


3/4 


A2 


0.277140 


0.120000 


0.105562 


0.044136 


0.823488 


A3 


0.170752 


0.035000 


0.402875 


0.089642 


0.897473 



Table 1: The four points on which the numerical results were collected for the anisotropic 
triangular lattice. The expected correlation exponent is also reported. 

above (see Appendix [B]). In four dimensions, the space orthogonal to a 2-D manifold is 
two-dimensional: thus we sampled the susceptibility along two directions w\ and W2 for each 
point, such that w\ _L W2 and Wi _L E. 

There is still an important caveat to keep in mind: while the above formulas provide the 
exact critical manifold, i. e. critical and integrable (satisfying the star-triangle relations), we 
know the complete critical manifold is here 3-dimensional. It is then to be expected that a 
particular direction in the plane spanned by (w\,W2) stays on criticality, but the chances of it 
being exactly one or the other of the two basis vector employed is vanishing; indeed, it turned 
out not to be the case for our numerical investigation. On the other hand, a direct trial-and- 
error numerical search for the orientation of the full critical manifold appears hardly feasible, 
time-consuming and of scarce conceptual relevance, so it has not been pursued further. 

The four points chosen for the numerical investigation, along with the corresponding 
prediction for correlation index, are listed in Table [TJ 

The (density of) energy susceptibility is defined as: 

C/ L (/3 1 ) = ((n-(n) /3i;L ) 2 )^ i , (5.35) 

where 

□ = ^2 ^2 {pi( a x<? y + T x r y ) + ai(a x a y T x T y )j , (5.36) 

that is, the average energy per link in the lattice according to the action for the system 
(the neighbouring pairs of sites (x,y)i are taken according to the corresponding anisotropic 
coupling). L is the side, in lattice sites, of the system, which has periodic boundary conditions 
and a square shape. Measurements are collected for different values of L and different points 
on the segments, and according to the predictions for a second-order transition the curves 
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for various system sizes should collapse one onto another if one plots 

F = L A -iu L (P) as a function of x = (ft - fflL$ , (5.37) 

with /3| denoting criticality (obviously, any other coupling would work). 

In this numerical investigation, 36 points were sampled along each of the two segments 
for the four points Al, AFZ, A2 and ^43; the collected statistics was of 36000 configurations 
per each data point, with six system sides ranging from L = 150 to L = 500. The algorithm 
used was an adaptation to the anisotropic Ashkin- Teller model of the standard cluster-based 
Swendsen-Wang prescription, based on identifying and flipping the Fortuin-Kasteleyn clusters 
alternatively in the a and r spin variables (with the other Ising spin playing the role of a 
local shift of the link coupling ft to ft ± on); for more details on the algorithm, we refer to 

The rescaled susceptibilities fall wery well to a universal curve F(x) in all four cases, 
thus confirming the location of the critical point along the segment and the accuracy of 
the prediction for the correlation index; to estimate the effect of systematic deviations on 
the latter (coming for instance from the finiteness of the system under study, or from the 
interplay between the anisotropy and the always square shape of the system), we compared 
the theoretical v with a value v* tuned as to yield the best collapse of the data, finding a 
discrepancy never exceeding about 5% (Fig. [6|). 

6 Conclusions 

In this paper we found the critical surface of the Ashkin- Teller model on the generic iso-radial 
graphs starting from the critical line of the anisotropic Ashkin- Teller model on the square 
lattice. Several special cases, such as the triangular and honeycomb lattice, were discussed 
in detail. Our method is in close connection with the duality relation and integrability of the 
critical surfaces. In other words, the critical surfaces that we considered in our work on the 
square lattice are self-dual and also integrable. Unfortunately they are not scanning the full 
critical surfaces in the anisotropic case: it seems, indeed, that there are regions of the full 
critical surface that are neither integrable nor self-dual even on the square lattice |37j . On 
the extracted critical surface for the Ashkin- Teller model on the generic iso-radial graph, we 
also found the free energy. 

In the special case of the anisotropic triangular lattice, we checked our conjecture nu- 
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Point A2, nu=0.80, segment W1 
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Figure 6: The susceptibilities for points Al, AFZ, A2 and A3 (top to bottom) along the two 
segments (left: w\; right: W2) as described in the text, rescaled according to Eq. I5.371 The 
plots employ the "experimental" correlation index v* as indicated in the figures' caption. 
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merically and found an excellent agreement with the conjectured behaviour. In addition we 
studied the lattice spinor variables for the Ising model in the non-critical case, finding the 
corresponding lattice equation of motions for Majorana fermions. It seems that a similar 
proceeding is possible for some particular points of the Ashkin- Teller model. Although we 
were not able to find the corresponding fermionic variables all along the critical line, we 
believe they exist. Another interesting study could be extracting the Thirring model, which 
is related to two Ising fermionic variables coupled by four-spinor interactions, with the same 
arguments; the difficulty in this direction might arise from writing down the multiplication 
of three disorder variables, which doesn't seem an immediate task. 

From the numerical point of view there are different ways to expand our study. For 
example, one can check the validity of our conjecture in the portion of parameter space 
beyond the Ising model, especially in the triangular-lattice case where it seems, in contrast 
with the square lattice, that one can reach asymptotically the X-Y point. Also, it could be 
interesting to explore the critical (anisotropic) Ising lines on different lattices; since they do 
not fall on the exactly solvable surfaces, there are no conjectures for their form, however it 
could be worth to see (at least numerically) their behaviour under a change of anisotropy 
angle. Another interesting direction is an investigation of the fractal properties of critical 
interfaces in the isotropic Ashkin- Teller model [381 EH] • 

A Ashkin-Teller model on triangular and honeycomb lattices 

In this Appendix further details about special points on the (isotropic) triangular- and 
honeycomb-lattice critical lines are provided. A startling feature of the phase space is that 
the boundaries of the physical region, i. e. that with non-negative Boltzmann weights, do not 
respect universality, so that different geometries allow to span different ranges of universality 
classes within the physical region (Fig. [5]). In the x\-X2 plane, this means that the physical 
threshold, X2 = —1 + 2xi, crosses the three critical lines, Eqs. 14.71 14.101 and 14.131 a t t ne al- 
ready mentioned points T n , XY^ and respectively, each one in its own universality class. 
In the f3, a plane, this means that those are the rightmost points that can be expressed with 
real values of the couplings (by sending them to infinity while keeping the sum (3 + a fixed 
to a specific value). As already observed (Eq. 12.20"]) . any attempt at translating non-physical 
points to this plane yields complex couplings. 

The most advantageous case is the triangular lattice: here, the whole critical line falls in 
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A point 


/3, a coordinates 


xi, X2 coordinates 


4-state Potts 


/3 = a = ^ 


X\=X2 = \ 


F-Z 


a _ Iog(2V3-3) _ log 3 
t ) 4 1 UL g 


Xl - i^jf-V f X2 _ 3 2V 2 


Ising 


a = , /3=M3 


Xi=2-y/3, X 2 = 7- Ay/3 


XY 


/3 = 12|2 _ a +OQ 


xi = y/2- 1 , x 2 = 2y/2- 3 



Table 2: Relevant critical points for the triangular isotropic lattice. 



O point 


0, a coordinates 


x\,X2 coordinates 


4-state Potts 


f3 = a = ^ 


Xi = x 2 = \ 


F-Z 


a _ log(3-2^2) _ log(2-v^) 
H — 4 ; u — 4 


xi = ^ , x 2 = V3(2 - x/3) 


Ising 


a - , f} - los(2 2 + ^ 


Xl = 73 ' X2 = 3 


U 


/3 - los(1 2 +v/5) a +oo 


xi = 4p±, x 2 = V5-2 


XY 




Xl = 2 ) ^2 = 



Table 3: Relevant critical points for the honeycomb isotropic lattice. 



the physical region, the XY point lying exactly on the boundaries; the coordinates for the 
relevant points are listed in Table [2j 

The square lattice has already been presented: it falls within the physical region only up 
to the T point, Eq. 12.171 The situation worsens for the honeycomb lattice: there, the critical 
line falls off the physical region starting from a point labelled U, that realises the following 
universality class: A^ = |-7r, u^j = | and r^j = y|. Table [3] contains the coordinates for the 
interesting points in this geometry. 

Duality provides correspondencies between points on the triangular and honeycomb lat- 
tices sharing the same universality class; on the other hand, the boundary of the physical 
region cuts the various critical lines in different points: this creates an interesting situation, 
in which a range of critical points (i. e. those with A > A^f) fall in the physical region on the 
triangular lattice but do not on the dual honeycomb realisation^- As already seen, the latter 
would correspond to values of /3, a with a nonzero imaginary component Nevertheless, 
duality must be respected even though, apparently, the Boltzmann weights are not always 

4 More precisely, this happens for all points in the (triangular-lattice) phase space, including the ones off 
criticality, that lie below the curve a = — i log cosh 2/3 but still within the physical region. 
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positive numbers. Indeed, a careful investigation reveals that the link-specific Boltzmann 
weights combine in such a way that, due to the geometric shape of the lattice and some 
constraints on the allowed spin configurations, the total weight for a given configuration is 
always a real number, and the resulting partition function is a well-defined, real and positive 
object (apart from a trivial overall prefactor). 

This situation can also be considered from a somewhat reversed point of view: starting 
from a statistical model on the honeycomb lattice with complex couplings and therefore ill- 
defined Boltzmann weights, by making use of duality one can re-express all of the observables 
in a triangular-lattice version of the model, with well-defined, real and positive Boltzmann 
weights; moreover, the latter version of the original problem would now be open to numerical 
investigation. That would be the case, for instance, if one wanted to study the point with 
r 2 = |, A = |7r, v = |, i. e. the one corresponding to the twisted M = 2 supersymmetric 
CFT; the triangular lattice is, moreover, the only way to perform numerical investigation in 
the vicinity of the XY point, which is accessible exclusively on this regularisation. 

B Anisotropic triangular lattice: complementary formulas 

This Appendix provides some more formulas specific for the anisotropic triangular lattice with 
two equal sides (Section [5]). From the parametrisation for the critical manifold, Eqs. 15.331 
and 15.341 two linearly independent vectors tangent to the surface in a given point P can be 
evaluated as: 



where the vectors are written as (/3i, a±, 02, o^)- The segments used to probe the susceptibility 
are then constructed along two orthogonal solutions wi,W2 to the system: 




(B.38) 



(B.39) 



ti ■ W = t2 ■ W = . 



(B.40) 



The above tangent vectors are explicitly given by: 




(B.41) 



(B.42) 
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where the following shorthands have been employed: 

d{8i) = cosh 2 2/3i; (B.43) 

, n . 1 — sinh 2 28i + 2 sinh 2 2a\ — 2 sinh 2ai cosh 2a\ . 

= m ; (B.44) 

m s~*t\ sinh 2G sinh 2«i 
«(A,ai,G) = sinh2/?i ; (B-45) 

r(ai) = 2(2 sinh 2ai cosh 2ai - 2 sinh 2 2ai - 1) . (B.46) 

Another complementary observation to this particular case of anisotropy concerns the 
limits 8\ — > and /?2 — >• 0, taken by staying on a fixed universality class; the former yields 
82 —> oo, while the latter just gives 8\ —> 8f , that is, the corresponding square-lattice 
couplings. This can be easily understood by considering that the two limits implement the 
removing of two or one of the sides of the triangle, respectively. In other words, the first limit 
gives a system that is just a bundle of disconnected one-dimensional systems (never getting 
a nonzero magnetisation), and the second results in a square lattice with isotropic couplings. 

We now provide three explicit "iso-A" lines within the critical manifold; the four-state 
Potts model is given by: 

ai = 8i , 8 2 = ^arctanh(l - 2tanh(2/3i)) ; (B.47) 



the F-Z line is expressed as: 



1 1 / 1 - sinh 20i J 2 + sinh 2 28i \ ,^ , N 

sinh 2a { = —= sinh 2$ , 82 = -arctanh v „ — ) ; (B.48) 

V2 2 V cosh 2B\ f 



and finally the Ising sub-manifold is: 



on = , 8 2 = ^arctanh f 1 Sll f 2 ^ ) . ( B .49) 
2 \ cosh 2a\ 



These are found by imposing A = 0, j, | in Eqs. 15.331 and 15.341 by means of Eq. 14.171 
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